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Abstract — Recent years have seen significant research in 
finding closed form expressions for the delay of the RLC 
interconnect which improves upon the Elmore delay 
model. However, several of these formulae assume a step 
excitation. But in practice, the input waveform does have 
a non zero time of flight. There are few works reported so 
far which do consider the ramp inputs but lacks in the 
explicit nature which could work for a wide range of 
possible input slews. Elmore delay has been widely used 
as an analytical estimate of interconnect delays in the 
performance-driven synthesis and layout of VLSI routing 
topologies. However, for typical RLC interconnections 
with ramp input, Elmore delay can deviate by up to 100% 
or more than SPICE computed delay since it is 
independent of rise time of the input ramp signal. We 
develop a novel analytical delay model based on the first 
and second moments of the interconnect transfer function 
when the input is a ramp signal with finite rise/fall time. 
Delay estimate using our first moment based analytical 
model is within 4% of SPICE-computed delay, and model 
based on first two moments is within 2.3% of SPICE, 
across a wide range of interconnects parameter values. 
Evaluation of our analytical model is several orders of 
magnitude faster than simulation using SPICE. We also 
discuss the possible extensions of our approach for 
estimation of source-sink delays for an arbitrary 
interconnects trees. 

Keywords- Delay Calculation, RLC Interconnect, 
Moment Matching, Ramp Input. 



I. 



Introduction 



'Accurate calculation of propagation delay in VLSI 
interconnects is critical for the design of high speed 
systems. Transmission lines effects now play an 
important role in determining interconnect delay and 
system performances. Existing techniques are based on 
either simulation or (closed-form) analytical formulas. 
Simulations methods such as SPICE give the most 
accurate insight into arbitrary interconnect structures, 
but are computationally expensive in total IC design 
stages. Faster methods based on moment matching 
techniques are proposed [8-9] but are still too expensive 
to be used during layout optimization. Thus, Elmore 
delay [2], a first order approximation of delay under 
step input, is still the most widely used delay model in 



the performance driven synthesis of clock distribution 
and Steiner global routing topologies. However, Elmore 
delay cannot be applied to estimate the delay for 
interconnect lines with ramp input source; this 
inaccuracy is harmful to current performance-driven 
routing methods which try to determine optimal 
interconnect segment lengths and widths (as well as 
driver sizes). Previous moment-based approaches can 
compute a response for interconnects under ramp input 
within a simulation-based methodology [10], but to the 
best of our knowledge there is no such analytical 
explicit delay estimation models proposed which is 
based on the first few moments. Recently, [3] presented 
lower and upper bounds for the ramp input response; 
their delay model is the same as the Elmore model for 
ramp input. Delay estimates for the analytical ramp 
input model are off by as much as 50% from SPICE- 
computed delays for 50% threshold voltage [10], and 
the analytical ramp input model cannot be used to 
obtain threshold delay for various threshold voltages. 
Ismail et. al. [4] used Elmore delay as an upper bound 
on the 50% threshold delay for RLC interconnection 
lines under arbitrary input waveforms. In this paper, we 
have proposed a novel and accurate analytical delay 
estimate for distributed RLC tree interconnects under 
arbitrary ramp input. To experimentally validate our 
analysis and delay formula, we model the interconnect 
lines having various combinations of sources, and load 
parameters, apply different input rise times, and 
obtained the delay estimates from SPICE, Elmore delay 
and the proposed analytical delay model. Over a wide 
range of test cases, Elmore delay estimates can vary by 
as much as 100% from SPICE computed delays. As the 
input rise time increases, Elmore delay deviates even 
further from SPICE results. 



II. 



Basic Theory 



A. Central Moments and Transmission Line Response 

For a simple input source terminated transmission 
line we can write the transfer function as, 



H (s )= 



V" _ , 

V (s) sC {R msk(fd)+Z sinh {yd))+{R 12 )swh(yd)+cosh[yd) 
i L s o so 



(1) 
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Where y = y(R + sL)(sC) is the propagation 
constant and Z =-\l{R + sL)l{sC) [ s the 
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characteristic impedance of the line. R, L and C are the 
per-unit-length resistance, inductance, and capacitance 
of the transmission line, respectively, d is the length of 
the line. The series resistance is given by, 
°i = "rf;-" 1 " °ter where R* is the driver resistance and 
R K r the termination resistance. We assume that at the 
given frequency of interest, the dielectric loss and 
conductance values are negligibly small. The driver 
resistance is assumed to be linear. Now the RLC 
interconnect can be considered as either lossless or 
lossy. 

B. Lossless Interconnect 

For an unloaded lossless transmission line driven by 
a step input, it is well known that the optimal 
termination resistance is R S =Z . With this termination, 
the ideal signal is the input step delayed by the time-of 
flight along the line, and is given by, T y=V LC a 
The following discussion shows that this ideal response 
is indeed obtained when the central moments of the 
impulse response are minimized. For the lossless line in 
Fig. 1, the transfer function is given by [5-6], 

1 

; R J Z J sinh (yd) + cosh (yd) (2) 

Where y and Z are the propagation constant and the 
characteristic impedance, respectively and are defined 
as, y=s^LC and Z =\ L/C . For this transfer 
function, the second and third central moments of the 
impulse response are symbolically given as: 

fi=-CLd 2 +R ,C 2 d 2 andfi= -2R C 2 Ld 3 +2R 3 C 3 d 3 
* s J s s 

(3) 



H(s) 



and localized (zero dispersion) about its mean, 
/u=4 LC d • Conversely, forcing the impulse 
response to be symmetric and localized about the mean 
ensures critical damping. 
So from (5) we can write the following equation: 



-sT , 



V (s)=V ( (s)e 
In case of ramp input, 

Vis)= V 



(6) 



DD 



(7) 



Substituting (7) in (6) we get, 
Vis) 



DD sT f 

e 



(8) 



Taking inverse Laplace transform of (8), 



v (t)=v DD (t- 



T f )u(t-T , 



(9) 



For calculation of the time delay we take V (s) = 0.5V D d 
at time t = T D and hence substituting in (9), we have, 



0-5V DD =V DD (T D - 



T f )u(t-T 



f> 



So for t > T f me x D is given as, 



T D =T f + 0.5 



(10) 



(11) 



The above equation (11) is our proposed closed 
form expression for delay for lossless transmission line 
RLC interconnects system. 

C. Lossless Interconnect Considering Mutual 
Inductance 

In order to calculate the exact time delay in two 
parallel RLC line, we consider the mutual inductance 
between two inductors as M. Fig. 2 shows two highly 
coupled transmission line system. 



"*(* 







Figure 1 . Lossless Transmission Line 



Solving for /'2 = ^ from (3) yields V L/C an d 

— V L/C as roots. Again solving /"3 = from (3) 

yields 0, 4 L/C and — V L/C as roots. The positive 
root provides the solution R s = Z , Then, the transfer 
function given as, 

H(s)- . . , - — — = e s f (4) 



sinh iyd) + cosh (yd 



Where T , 



j — V LC d j s me time of flight. Then it 

is easy to show that this transfer function provides the 
desired ideal waveform at the output of the transmission 
line 

v o (t) = Vi (t-T f ) (5 ) 

From (5), it can be inferred that the ideal impulse 
response for a lossless transmission line is symmetric 




victim driver 

r 

aggressor driver 

rVvV 

V *(T) "* 



Figure 2. Highly Coupled Transmission Line 

1) Mutual Inductance 
The mutual inductance M of two coupled 
inductances Li and L 2 is equal to the mutually induced 
voltage in one inductance divided by the rate of change 
of current in the other inductance 



M=EJ 



di j 
, dt 
(di, 



(12) 



(13) 



If the self induced voltages of the 
inductances Li and L 2 are Ei, and E 2s , respectively, for 
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the same rates of change of the current that produced 
the mutually induced voltages Ei m and E 2m , then: 



M= 


(e ? \ 

2m 




fc'i 




\ Is/ 




*1 \ 


M- 


lm 




\ 2s I 



*L X 



*L, 



Combining equations (14) and (15) we get, 



M- 



E\ m E 2m 

E \s E 2s 



4lTl : 



M 



L\ L 2 



1/2 



(14) 



(15) 



(16) 



Where k M is the mutual coupling coefficient of the 
two inductances Li and L2. If the coupling between the 
two inductances Li and L 2 is perfect, then the mutual 
inductance M is: 
M = (UL 2 y/2 
2) Odd Mode 

When two coupled transmission lines are driven 
with voltages of equal magnitude and 180° out of 
phase with each other, odd mode propagation occurs. 
The effective capacitance of the transmission line will 
increase by twice the mutual capacitance, and the 
equivalent inductance will decrease by the mutual 
inductance [7]. In Fig. 3, a typical transmission line 
model is considered where the mutual inductance 
between aggressor and victim connector is represented 
as M12. Li and L 2 represent the self inductances of 
aggressor and victim nodes, respectively, while C c , C 
denotes the coupling capacitance between aggressor and 
victim and self capacitance, respectively. 



3) Even Mode 
When two coupled transmission lines are driven 
with voltages of equal magnitude and in phase with 
each other, even mode propagation occurs. In this case, 
the effective capacitance of the transmission line will be 
decreased by the mutual capacitance and the equivalent 
inductance will increase by the mutual inductance [7]. 
Thus, in even-mode propagation, the currents will be of 
equal magnitude and flow in the same direction. The 
magnetic field pattern of the two conductors in even- 
mode is shown in Fig. 5. The effective inductance due 
to even mode of propagation is then given by, 



even 1 2 



(18) 




Figure 5. Magnetic Field in Even Mode 

III. Proposed Delay Model 

A. Calculation of Delay for Even Mode 

From (1) and for a simple input source terminated 
transmission line, we can write the transfer function as, 



His 



' r : ( s ) 
l 



sCjR s oo^ly e d) + Z tV imhly 1 ,d})+{R s !Z x )imh{y t ,d)+<:oi\i(y r d) 




Figure 3. An Example for Two Parallel Transmission Line Model 

Assuming that Li = L 2 = L , the currents will be of 
equal magnitude but flow in opposite direction [7]. 
Thus, the effective inductance due to odd mode of 
propagation is given by, 

A*/</ = A ~ ^2 (17) 

The magnetic field pattern of the two conductors in 
odd-mode is shown in Fig. 4. 




Li, 

(19) 

Where is the y e =-J(R + s{L+ M)){sC) i s the 

propagation constant for even mode and 

Z oe =-J(R + s(L+M))/(sC) i s the characteristic 

impedance for even mode of the line. R, L, M and C are 
the per-unit-length resistance, inductance, mutual 
inductance and capacitance parameters of the 
transmission line, respectively, d is the length of the 
line, and the series resistance is given by 

™s = dr~^~^ter where Rdris the driver resistance and 
R tcr the termination resistance. We assume that the 
dielectric loss and hence the conductance, G to be 
negligibly small. The driver resistance is assumed to be 
linear. For an unloaded lossless transmission line driven 
by a step input, it is well known that the optimal 
termination resistance is R s =Z oc . With this termination, 
the ideal signal is the input step delayed by the time-of 
flight along the line, is given by, 

T ,=^{L + M)Cd _ 



Figure 4. Magnetic Field in Odd Mode 



The following discussion shows that this ideal 
response is indeed obtained when the central moments 
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of the impulse response are minimized. For the lossless 
line in Fig.l, the transfer function is given by, 



H s)=- 



1 



R IZ )sinh(v <i) + cosh(y d) 
s oe ' e w e ' 

(20) 
Where y e =syJ(L+M)C an( j 

Z oe —y(L + M)/C for this transfer function, the 

second and third central moments of the impulse 
response are symbolically given as: 

M 2 = -C(L+M)d 2 +R 2 C 2 d 2 

S 

f.i 3 = -2R^C 2 (L + M)d 3 +2R 3 C 3 d 3 

s 

oil 



6 6 



(21) 

Solving for M2~® from equation (21) yields 

^{L + M)/C and -V(Z+M)/C as roots. Again 

solving /"3 = from equation (21) yields 0, 

^(L + M)/C and -AL+M)IC as roots The 

positive root provides the solution R s = Z„ e , Then, the 
transfer function given as, 



H(s) = 



- S T 



sinh(y d)+cosh(y d 



-—e 



(22) 

Where, T f =^l{L+M)Cd i s the time of flight. 

Then it can be easily shown that this transfer function 
provides the desired ideal waveform at the output of the 
transmission line is: v o\ t )~ V A t ~ 1 fj 

From above, it can be inferred that the ideal impulse 
response for a lossless transmission line is symmetric 
and localized (zero dispersion) about its mean, 
fi=y(L+ M)C d ■ Conversely, forcing the impulse 
response to be symmetric and localized about the mean 
ensures critical damping. 

So from equation (22) we can write the following 
equation: 

V (s)=V i (s)e ST '' (23) 

In case of ramp input, 
V, 



V : s =- 



DO 



Substituting (24) in (23) we get, 

77 I \ V DD ~ ST f e 
S 

Taking inverse lapalce transform of (25), 
V (t) = V DD (t-T fe )u(t-T fe ) 



(24) 



(25) 



(26) 



For the calculation of the time delay we take V (s) = 
0.5V D d at time r=T D and hence substituting in (26), we 
have, 



Q.5W DD =V DD {T D -T fe )u(t-T fe ) (27 ) 

So for t > T f e ( j D ^ given as, 
T D =T fe + 0.5 (28 ) 

The above equation (28) is our proposed closed 
form expression for delay for lossless transmission line 
RLC tree circuit in even mode and with mutual 
inductance. 
B. Calculation of the Delay in Odd Mode 

Again from equation (1) for a simple input source 
terminated transmission line we can write the transfer 

function as, 

r l') 



H(s)-- 



VAs) 



1 



S C l (R s cmh(yJ)+Z mS mii{ ro d))-t-{RjZj S Mi(y d)+co S h(y d) 



(29) 
Where, y =^l(R + s(L-M))(sC) ls the 



propagation 



constant 



and 

Z 00 = ^(R + s(L-M))/(sC) i s the characteristic 

impedance for odd mode of the line, respectively. R, L, 

M and C are the per-unit-length resistance, inductance, 

mutual inductance and capacitance parameters of the 

transmission line, respectively, d is the length of the 

R = R 4- R 
line, and the series resistance is given by s Jr "'' 

. Where R* is the driver resistance and R tcr the 

termination resistance. We assume the dielectric loss 

and hence the shunt conductance, G to be negligibly 

small. The driver resistance is assumed to be linear. 

For an unloaded lossless transmission line driven by 

a step input, it is well known that the optimal 

termination resistance is R s =Z o With this termination, 

the ideal signal is the input step delayed by the time-of 

flight along th e line, is given by, 

T fb =yJ(L-M )Cd 

The following discussion shows that this ideal 
response is indeed obtained when the central moments 
of the impulse response are minimized. For the lossless 
line in Fig. 1 , the transfer function is given by, 

1 



H(s)-- 



( R s lz oo)s^y d ) 



-cosh(y d) 



(30) 
Where y =S^j{L — M)C is the is the 

propagation constant and Z 00 =y(L — M)/C is the 

characteristic impedance for this transfer function, the 
second and third central moments of the impulse 
response are symbolically given as: 

fi 2 = -C(L-M)d 2 +R jC 2 d 2 



u =-2R C 2 (Z-MW 3 +2R 3 C 3 ^ 3 

r 3 s v ' s 



(31) 
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Solving for jJ. 2 — from equation (31) yields 

^{L-M)/C and -V(Z-M)/C as roots. 

Again solving jU^ = from equation (31) yields 0, 

^{L-M)IC and -V(I-M)/C as roots. The 
positive root provides the solution R s = Z 0o , Then, the 
transfer function may be expressed as, 



His) 



1 



-sT, 



(32) 



sinh(y o d)+cosh(7 o d) 

Where T f =4{L-M)Cd is the time-of- 

flight. Then it is easy to show that this transfer function 
provides the desired ideal waveform at the output of the 
transmission line is v {t) = v i [t — 1 ^ ) . 

From above, it can be inferred that the ideal impulse 
response for a lossless transmission line is symmetric 
and localized (zero dispersion) about its mean, 
/u=y(L — M)Cd conversely, forcing the impulse 
response to be symmetric and localized about the mean 
ensures critical damping. 

So from (32), we can write the following equation: 

V (s)=V t (s)e-"'- (33) 

In case of ramp input, 
VAs)=^ (34) 



Substituting (34) in (33) we get, 



V (s) 



V 



DD 



-sT, 



(35) 



Taking inverse Laplace transform of equation (35) 

V (t)=V DD (t-T fo )u(t-T fo ) (36) 

In order to calculate the time delay we take Vo(s) = 
0.5Vdd at time t = T D and hence putting in equation 
(36), we have, 



0.5V DD =V DD (T D -T fo )u(t-T fo ) 

So for t > T , Q the T D is given as, 



D 



-T fo + 0.5 



(37) 



(38) 



The above equation (38) is our proposed closed 
form expression for delay for lossless transmission line 
RLC tree circuit in Odd mode and with mutual 
inductance. 

IV. Experimental Results 

In the case of very high frequencies as in GHz scale, 
inductive effect comes into the important role and it 
should be included for complete delay analysis. The 
configuration of circuit for simulation is shown in 
Figure 6. 



TtoiffliiiMiliiie Jrarntnmionte 

I 1 BU3 M B1C 1 1 

to TiaiMiMMiline TrarnttiMoii line Tianminim line 

-^s~- 1 rh I — I 110 — I B 



~t 



TrarutninionliiK Tmimiiiioii line 

- I S H S i—i 



Figure 6. An RLC Tree Example 

For each RLCG network source we put a driver, 
where the driver is a step voltage source followed by a 
resistor. The results are based on equation (38) for 0.18 
um process. The left end of the first line of Fig. 6 is 
excited by IV ramp form voltage with rise/fall times 0.5 
ns and a pulse width of Ins. In table 1, the 50% delay 
for even mode and the Elmore delay is compared for 
various values of the driver resistance Rd and the load 
capacitance Cl when the length of the RLC interconnect 
is kept constant. In the similar way, in table 2 the 50 % 
delay for odd mode and the Elmore delay are compared. 

Table I 
Experimental result under Ramp input for even mode 



Ex 


Rs 
(£}) 


Cl(iT) 


L(um) 


Ted 
(ps) 


Proposed Delay 
Model (ps) 


1 


1 


10 


100 


0.1251 


0.1342 


2 


2 


50 


100 


0.1567 


0.1576 


3 


5 


750 


100 


0.4589 


0.4765 


4 


10 


1000 


100 


0.9310 


0.9142 


5 


50 


1500 


100 


0.3920 


0.3675 


6 


100 


1500 


100 


0.5955 


0.5762 


Table II 
Experimental result under Ramp input for odd mode 


Ex. 


R s 
(£2) 


Cl(iT) 


L(iun) 


Ted 
(ps) 


Proposed Delay 
Model (ps) 


1 


1 


10 


100 


0.1065 


0.1127 


2 


2 


50 


100 


0.2597 


0.2376 


3 


5 


750 


100 


0.4943 


0.4869 


4 


10 


1000 


100 


0.9876 


0.9792 


5 


50 


1500 


100 


0.3724 


0.3684 


6 


100 


1500 


100 


0.5732 


0.5989 



In table 3 and table 4, comparative result of our 
proposed model delay with the SPICE delay are given 
in the similar way as we did for the comparison of our 
proposed model and Elmore delay model discussed 
above as in table 1 and table2. 

V. Conclusions 

In this paper we have proposed an accurate delay 
analysis approach for distributed RLC interconnect line 
under ramp input. The use of transmission line model in 
our study gives a very accurate estimate of the actual 
delay. We derived the transient response in time domain 
function of ramp input. We can see that when 
inductance is taken into consideration, the Elmore 
approach could result an error of average 10% 
compared to the actual 50% delay calculated using our 
approach. 
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Table III 
Comparative result with spice under Ramp input for even mode 



R, 
(O) 


C L (fF) 


L(um 

) 


SPICE 

(PS) 


Proposed Delay 
Model (ps) 


% 
Error 


1 


10 


100 


0.1451 


0.1342 


7.51 


2 


50 


100 


0.1595 


0.1576 


1.19 


5 


750 


100 


0.4789 


0.4765 


0.51 


10 


1000 


100 


0.9317 


0.9142 


1.87 


50 


1500 


100 


0.3928 


0.3675 


6.44 


100 


1500 


100 


0.5997 


0.58 


3.91 



Table IV. 

Comparative result with spice under Ramp input for odd mode 



Rs 

(Q) 


CL(f 
F) 


L(um 

) 


SPICE 
(PS) 


Proposed Delay 
Model (ps) 


% 
Error 


1 


10 


100 


0.1185 


0.1127 


4.89 


2 


50 


100 


0.2588 


0.2376 


8.19 


5 


750 


100 


0.4994 


0.4869 


2.5 


10 


1000 


100 


0.9843 


0.9792 


0.51 


50 


1500 


100 


0.3792 


0.3684 


2.84 


10 



1500 


100 


0.5787 


0.5989 


3.49 
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